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I.  INTRODUCTION 


The  general  scope  of  this  project  is  to  investigate  properties  of  ferroelectric  ceramics 
using  new  measurement  techniques  and  to  look  for  new  technological  uses  of  these  materials  in 
device  applications.  One  main  focus  of  the  effort  is  a  study  of  interfaces  between  counterpoled 
regions  of  a  ferroelectric  ceramic.  Such  counterpoled  ceramics  have  been  used  in  the  past 
for  focused  Fresnel  zone  plate  transducers,  bimorphs  and  tuning  fork  resonators.  In  such 
applications  it  is  important  to  know  the  structure  of  the  interface  and  how  it  relates  to 
material  composition  and  history. 

Powerful  techniques  for  experimental  study  of  these  questions  involve:  (I)  elastic  wave 
scattering  by  an  interface,  (2)  elastic  wave  guidance  along  an  interface,  (3)  laser  probe 
pyroelectric  measurements  of  the  electric  polarization  profile  through  an  interface,  (4)  laser 
probe  measurements  of  the  mechanical  vibration  profile  through  an  interface.  Items  (1)  and 
(2)  lead  to  potential  device  applications  in  the  form  of  grating  filters  and  interface  delay 
lines.  A  theoretical  basis  for  treating  uniform  arrays  of  interfaces,  such  as  might  be  used 
in  these  devices,  requires  a  general  analysis  of  high  frequency  elastic  wave  propagation  in 
spatially  periodic  media.  Such  a  theory  can  be  constructed  by  using  space  harmonic  analysis, 
as  developed  for  the  vibration  spectra  of  crystal  lattices,  and  is  also  applicable  to  periodicities 
in  one-,  two-  or  three-dimensional  dimensions.  These  multidimensional  periodicities  occur  in 
a  variety  of  PZT-epoxv  composite  materials,  and  a  theory  of  elastic  wave  propagation  is  of 
fundamental  importance  in  understanding  this  class  of  ferroelectric  materials. 

The  material  interfaces  discussed  above  may  be  introduced  extrinsically  by  inserting 
actual  material  inclusions  in  a  matrix  medium  (as  in  the  case  of  the  family  of  PZT-epoxv 
composite  transducer  materials),  or  intrinsically  by  altering  the  material  properties  with 
selectively  applied  spatial  distributions  of  electric  field  bias.  In  the  second  case,  one  can 
distinguish  between:  (1)  intrinsic  material  inhomogeneities  introduced  by  poling  patterns 
of  electrical  polarization  into  remanent-type  ferroelectrics  (such  as  PZT  or  PLZT),  and  (2) 
intrinsic  inhomongeneities  introduced  by  selective  electric  field  biasing  of  relaxor-type  fer¬ 
roelectrics  (such  as  PMN).  Materials  with  intrinsic  inhomogeneities  offer  exciting  possibilities 
for  realizing  a  variety  of  planar  signal  routing  and  processing  devices  on  small  ferroelectric 
wafers,  fabricated  by  simple  inexpensive  poling  and  biasing  techniques.  Such  devices  may 
be  either  slowly  or  rapidly  switched  and  programmed,  according  to  the  choice  of  either  a 
remanent-type  or  a  relaxor-type  ferroelectric  material. 
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n.  COMPOSITE  RESONATOR  AND  TRANSDUCER 

STUDIES 

These  studies  are  being  pursued  in  collaboration  with  R.T.  Gururaja,  W.A.  Schulze  and 
L.E.  Cross  at  the  Materials  Science  Laboratory  of  the  Pennsylvania  State  University,  and 
with  Dr.  W.  Smith  at  the  Philips  Research  Laboratory.  The  initial  goal  was  to  investigate  the 
high  frequency  dynamic  behavior  of  3:1  connected  PZT-epoxy  composites  in  plate  transducer 
configurations,  to  better  understand  the  physics  of  these  materials,  and  to  explore  their 
exploitation  for  applications  other  than  hydrophones.  Progress  made  in  formulating  a  general 
Floquet  (or  space  harmonic  expansion)  theory  of  elastic  vibrations  in  composites  now  makes 
it  possible  to  extend  the  goals  of  this  research.  The  dynamic  behavior  of  composite  materials 
consisting  of  general  periodic  lattices  of  inclusions  in  a  matrix  can  be  analysed  using  this 
theory,  which  has  been  tested  in  certain  specific  configurations  against  experimental  results.  A 
capability  now  exists  for  quantitatively  studying  and  exploiting  periodic  composite  materials 
(whether  piezoelectric  or  nonpiezoelectric)  in  meeting  Navy  needs  for  improved,  more  reliable 
and  cost-effective  techniques  for  detecting,  generating,  and  controlling  acoustic  vibrations 
in  liquids  and  solid  structures — including  specifically  the  development  of  anechoic  structures 
and  the  control  of  undesirable  noise  in  detectors. 

A.  Theory 

Basic  features  of  the  theory  were  outlined  in  the  Annual  Progress  Report  for  1  February 
1982-31  January  1983,  and  initial  encouraging  comparisons  with  experiment  were  reported 
there.  During  the  year  covered  by  the  present  report,  these  results  were  refined  and  compared 
in  more  detail  with  electrical  transducer  impedance  measurements,  time  delay  measurements 
of  acoustic  velocity,  and  laser  probe  vibration  pattern  measurements.1-3  References  l  to  3 
described  the  theory  for  unbounded  isotropic  composites  with  two-dimensional  periodicity 
in  the  xy  plane  (summarized  in  Appendix  A).  Explicit  calculations  have  been  made  for  the 
■‘superlattice'’  unit  cell  geometries  shown  in  Fig.  1. 

Although  the  general  formulation  is  capable  of  treating  Floquet  waves  propagating  in 
arbitrary  directions,  attention  was  initially  focused  on  wave  propagation  either  normal  to  the 
plane  of  the  periodicity  ( xy-plane)  or  parallel  to  that  plane.  These  propagation  directions  cor¬ 
respond.  in  the  first  case,  to  the  longitudinal  thickness  resonances  of  a  composite  transducer 
and.  in  the  second  case,  to  lateral  resonances  of  the  transducer.  Lateral  resonances  occur 
in  all  transducers.4  However,  certain  resonances  in  composite  transducers  are  enhanced  by 
Bragg-scattering  from  the  lateral  periodicity,  and  by  phase-matched  coupling  of  the  electri¬ 
cal  excitation  to  selected  lateral  resonances  through  the  periodically  distributed  piezoelectric 
coupling. 

Figure  2(a)  gives  a  typical  input  impedance  versus  frequency  curve,  measured  on  a  Penn 
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FIGURE  1 

Two-dimensional  unit  cell  geometries.  Area  of  unit  cell  is  A.  Area  of  PZT  (shaded)  is  a. 


State  sample  and  showing  three  strong  resonances.  Here,  ft  is  the  longitudinal  thickness 
resonance  and  fa.  /£2  are  the  two  dominant  lateral  standing  wave  resonances — along  the  unit 
cell  edge  and  the  unit  cell  diagonal,  respectively  (Fig.  2(b)).  Since  the  sample  is  uniformly 
electroded,  the  piezoelectric  elements  (shown  stippled  in  Figs.  1  and  2(b))  are  driven  in  phase, 
so  that  the  elastic  vibrations  excited  must  also  be  in  phase  at  these  points.  This  means  that 
the  lateral  resonances  observed  in  Fig.  2(a)  consist  of  standing  waves  with  one  full-wavelength 
spacing  between  the  piezoelectric  elements — along  the  unit  cell  edge  for  fa  and  along  the 
unit  cell  diagonal  for  /£2.  Such  standing  wave  patterns  are  shown  schematically  in  Fig.  3  ((a) 
for  fa  and  (b)  for  /(2).  These  patterns  will  be  discussed  more  fully  in  Section  II. B  below,  in 
connection  with  laser  probe  measurements  of  the  motion  of  composite  resonators. 

Our  study  of  lateral  resonances  in  composite  transducers  was  initially  approached  by 
assuming  wave  propagation  in  the  ry-plane  of  an  infinite  composite.  The  model  was  further 
simplified  by  supposing  that  only  a  r-polarized  elastic  displacement  existed,  reducing  the 
elastic  wave  function  to  a  scalar  form  similar  to  that  appearing  in  electronic  band  theory.  In  a 
first  approximation,  coupled  mode  theory  (which  retains  only  simultaneously  ■‘resonant”  space 
harmonics)  is  used  to  find  the  edges  of  the  stop  bands  for  wave  propagation  in  the  composite 
medium.  Stopband  edges  are  located  at  frequencies  where  Bragg-scattering  resonances  occur 
between  certain  planes  of  PZT  rods  (Bragg  planes)  in  the  periodic  array  of  rods.  (Resonances 
occur  at  frequencies  for  which  the  Bragg  planes  are  spaced  by  integral  multiples  of  X/‘2  for 
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FIGURE  2 

Resonance  spectrum  of  a  periodic  composite  plate  resonator:  (a)  Typical  measured  impedance 
curve:  (h)  Corresponding  standing  wave  patterns. 
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FIGURE  3 


Unit  cell  geometries  and  laser  scan  paths  for  d  =  X  Bragg  resonances:  (a)  Standin 
along  cell  edge  |/(1 ):  (b)  Standing  waves  along  cell  diagonal  (ft 2). 


the  2-polarized  shear  wave  assumed  in  the  model.)  These  “Bragg”  resonances  are  just  the 
enhanced  lateral  resonances  referred  to  above. 

At  X/2  spacing  between  the  Bragg  planes,  the  resonant  wave  functions  are  standing  waves 
along  the  unit  cell  edge,  the  unit  cell  diagonal,  and  perpendicular  to  the  higher  Bragg  planes. 
These  resonances  cannot,  however,  be  excited  piezoelectrically  in  a  uniformly  electroded 
sample  because  neighboring  PZT  rods  vibrate  180°  out  of  phase.  The  lowest  lateral  mode 
that  can  be  excited  piezoelectrically  in  this  geometry  is  that  which  occurs  when  the  edge  of 
the  unit  cell  equals  one  full  wavelength  (fti  in  Fig.  2).  At  this  resonant  frequency  several 
allowed  Bragg-scattering  resonances  occur,  but  only  the  mode  shown  in  Fig.  3  has  the  PZT 
rods  all  vibrating  in  phase.3  For  the  longitudinal  thickness  resonance  (fe  in  Fig.  2),  Bragg- 
scattering  resonances  do  not  appear,  because  there  is  no  periodicity  of  the  medium  along  the 
r-propagation  direction.  In  this  case  the  simplest  model  is  a  2-propagation  longitudinal  wave 
with  properties  determined  by  p  and  ?33  (spatial  averages  over  the  unit  cell),  and  with  weakly 
excited  space  harmonics  because  no  resonance  of  the  space  harmonics  can  occur.  In  other 
words,  the  Bragg-scattering  condition  cannot  be  satisfied. 

Comparison  of  experimentally  measured  resonant  frequencies  with  those  calculated  from 
the  simple  model  described  above  gives  good  agreement  for  the  longitudinal  thickness  mode 
but  only  fair  agreement  for  the  lateral  modes.1.  One  reason  for  this  is  that  the  theory 
uses  a  spatial  average  of  the  density  and  the  elastic  stiffness  constants  cu  (see  Eqs.  (2)-(4) 
in  Appendix  A)  to  model  the  effective  properties  of  the  medium  at  low  frequencies,  where 
scattering  from  the  periodicity  is  unimportant.  This  corresponds  to  the  constant  strain  model 
for  the  static  (or  low  frequency)  behavior  of  a  composite.  If  a  particular  strain  component 
(the  7-th)  has  the  same  value  (5y)  in  the  plastic  matrix  and  in  the  PZT  rod,  the  average  /-th 
component  of  stress  over  the  unit  cell  is  (see  Fig.  1) 


^(-4  —  g){cij)piattie  +  fl(c/y)pzr^^  _  ■^lJgJ 


(1) 


corresponding  to  the  spatial  averaging  of  the  stiffness  c//  in  Eq.  3(a)  of  Appendix  A.  For  the 
longitudinal  thickness  resonance  (/<-)  the  plastic  matrix  and  the  PZT  rods  vibrate  like  springs 
in  parallel,  and  the  model  of  Appendix  A  reduces  to  this  constant  strain  description  at  low 
frequencies. 

For  the  lateral  resonances  in  Fig.  2(b),  the  plastic  and  the  PZT  are  elastic  members 
connected  in  series,  so  that  the  stress  (not  the  strain)  is  common  to  both;  and  the  low 
frequency  behavior  of  the  medium  is  best  represented  by  the  constant  stress  model.  In  this 
case  the  counterpart  of  Eq.  (1)  is 


^(A  -  a)(sij)piattie  +  g(^/y)pzr^j^  — 
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for  the  average  strain  over  the  unit  cell.  Vibrations  of  the  composite  material  at  low 
frequencies  are  now  most  accurately  calculated  from  the  spatially  averaged  compliances  su , 
rather  than  the  c/j  of  Appendix  A. 

Another  reason  for  the  observed  discrepancies  between  experiments  and  the  theory  of 
Appendix  A  is  that  the  stress-free  boundary  conditions  at  the  surfaces  of  the  plate  have  not 
been  considered.  This  point  was  discussed  briefly  in  Reference  2.  It  was  noted  there  that 
additional  coupling  among  the  space  harmonics  occurs  at  the  boundaries,  and  that  new  plane 
wave  terms  (not  contained  in  the  original  set  of  space  harmonics)  are  also  generated  by  elastic 
mode  boundary  coupling. 

Appendix  A  (excerpted  from  Reference  3)  summarizes  the  Floquet  theory  of  wave 
propagation  in  an  infinite  isotropic  elastic  composite  with  two-dimensional  periodicity.  This 
theory  is  not  adequate  for  a  complete  treatment  of  the  composite  plate  transducer  prob¬ 
lem,  because  it  does  not  include  plate  boundary  conditions.  As  was  noted  in  the  previous 
paragraph,  the  addition  of  boundary  conditions  to  the  theory  of  Appendix  A  introduces  a 
number  of  complications.  Several  attempts  have  been  made  to  follow  this  approach.  However, 
it  was  found  that  (because  of  additional  boundary  couplings  introduced  between  space  har¬ 
monic  orders  and  between  different  polarizations  of  the  same  space  harmonic  order)  this  is 
not  a  useful  wav  to  implement  the  coupled  mode  approximation  for  composite  plates  or  to 
iterate  the  coupled  approximation  up  to  solutions  of  increased  accuracy. 

A  solution  to  the  above  difficulty  is  found  by  taking  the  basis  functions  of  the  Floquet 
solution  to  be  Lamb  and  SH  waves  of  the  plate,  rather  than  plane  waves.  This  procedure 
has  been  carried  out  in  detail  for  the  case  of  a  composite  plate  with  material  parameter 
periodicity  along  one  direction  (x)  only  (Fig.  4).  Details  of  the  analysis  are  given  in  Appendix 
B.  where  the  general  set  of  coupled  space  harmonic  amplitude  equations  is  shown  to  be  of 
the  form 

(cJ  =  ^  ]  fc-mn.ni/Q'ni/  (3) 

ni s 

In  Eq.  (3)  rn,  n  are  one-dimensional  space  harmonic  indices,  and  p,i/  are  Lamb  mode  indices. 
(For  a  one-dimensional  periodicity,  the  Lamb  and  SH  wave  solutions  are  separable  by  sym¬ 
metry.  The  case  of  two-dimensional  periodicity  can  also  be  treated  by  the  same  method). 
Using  Eq.  (3).  the  coupled-mode  approximation  is  easily  implemented  by  selecting  only  the 
near  ‘‘resonance"  terms  on  the  left-hand  side.  However,  it  must  be  noted  that,  in  these  equa¬ 
tions.  simultaneous  "resonance"  (Bragg-scattering)  may  occur  either  between  space  harmonics 
of  the  same  Lamb  wave  species,  or  between  space  harmonics  of  different  Lamb  wave  species. 
This  adds  a  great  many  stop  bands  in  the  wave  dispersion  curves.  Nevertheless,  the  Lamb 
wave  approach  of  Appendix  B.  as  distinct  from  the  plane  wave  approach  of  Appendix  A,  does 
give  a  clear-cut  method  for  describing  and  studying  the  stop  band  (or  resonance)  structure 
of  the  composite  plate. 
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An  additional  feature  of  the  Lamb  wave  formulation  is  that  it  leads  to  the  average 
eomplian.ee  (constant  stress)  material  description  of  Eq.  (2).  rather  than  the  average  stiffness 
(constant  strain)  description  of  Eq.  (1).  It  therefore  more  accurately  represents  the  lateral 
resonances  in  a  composite  plate.  This  point  is  confirmed  by  the  comparison  of  theory  and 
experiment  in  subsection  II. B.  below 

B.  Experiment 

Comparisons  of  theory  and  experiment  have  been  made  using  resonance  spectra  of  the 
eleriric  input  impedance  (Fig.  2)  and  laser  probe  measurements  of  elastic  resonance  vibration 
patterns  o\er  a  unit  cell  of  the  composite.  The  second  class  of  measurements  has  proved  to  be 
a  powerful  tool  for  positively  identifying  and  classifying  resonances.  For  reasons  mentioned 
in  subsection  II. A.  only  the  resonances  f(.  Jn  and  jt2  of  Fig.  2b  are  strongly  excited  in  a 
uniformly  electroded  sample,  and  comparisons  are  made  for  only  these  resonances. 

Many  laser  probe  measurements  have  been  performed,  but  only  one  example  will  be  given 
here  (Fig.  5).  The  figure  shows  recorded  amplitude  and  phase  of  the  vibration  amplitude 
measured  along  a  path  between  the  rods  in  Fig.  3.  Experimentally  observed  amplitude 
maxima  in  Fig.  3  are  seen  to  occur  at  the  zero  position  and  at  point  6  of  Fig.  3.  In  Fig.  3(a). 
superposition  of  the  x-  and  y-directed  standing  waves  predicts  a  zero  at  the  zero  position 
and  a  maximum  at  point  6,  In  Fig.  3(b).  superposition  of  the  two  standing  waves  predicts 
maximum  amplitudes  at  both  the  zero  position  and  at  point  6.  The  latter  behavior  is  in 
agreement  with  the  experimental  result  in  Fig.  5.  It  can  also  be  seen  by  superposing  the 
two  standing  waves  in  Fig.  3(b)  that  the  vibration  maximum  at  the  zero  position  and  that 
at  point  b  arc  18(T  out  of  phase,  again  in  agreement  with  the  experiment  of  Fig  5.  This 
comparison  clearly  identifies  the  motion  as  belonging  to  the  diagonal  resonance  of  part  (b) 
m  Fig  3.  Corroboration  is  provided  by  measurements  taken  at  points  a.  c  and  d  in  Fig  3. 
which  were  also  found  to  be  consistent  with  the  superposed  standing  wave  patterns  in  Fig. 
3(  b  | 

The  above  laser  probe  technique  was  used  at  Sanford  to  identify  the  vibration  patterns 
nf  tht’  resonances  observed  in  certain  selected  samples  among  those  listed  in  Table  I.  For 
these  particular  samples,  the  resonance  spectra  (as  in  Fig.  2(a))  measured  at  both  Stanford 
and  Penn  State  were  found  to  he  in  close  agreement.  Most,  of  the  sample  frequencies  listed  in 
Table  I  were  measured  at  Penn  State  only.1  and  the  resonant  mode  identification  was  made 
by  comparison  with  samples  that  were  laser-probed. 

In  Table  [  the  piezoelectric  rods  are  of  PZT-5A  and  the  polymer  of  Spurrs  epoxy .  Material 
parameters  used  in  comparing  theory  and  experiment  are  p  =  7.75  X  103  Kg/m3,  Yt  = 
1032  m/s.  1,  =  2263  m/s  for  PZT-5A.  and  p  =  1  14  X  103  Kg/m3.  Y(  =  2207  m/s. 

I  ,  =  1003  m/s  for  Spurrs  epoxy.  The  parameters  for  PZT-5A  were  obtained  from  the 
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PHASE  (DEGREE)  AMPLITUDE  (dB) 


Laser  scan  measurement  of  acoustic  displacement  at  /t2,  scanned  between  the  rods  (Fig.  3(b)). 
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TABLE  I 


> 


Resonant  Modes  for  PZT  Rod-Polymer  Composites' 

(The  numbers  on  the  left  correspond  to  the  points  plotted  in  Figs.  6-11) 


Label  of 
Plot  ted 
Point 

Sample 

Number 

RJ  PZT 

Thickness 
b  (mm) 

Spacing 
d  (mm) 

fi  (anti- 

resonance) 

(KHZ) 

fz  (reson¬ 
ance) 
(KHZ) 

•  /3  (reson¬ 
ance) 
(KHZ) 

Frequency 

Plotted 

1 

128" 

10R? 

1.9 

1.27 

- 

804 

1096 

fz 

a 

a 

5°o 

2.59 

1.58 

- 

604 

786 

fz 

3 

128" 

10°F 

1.9 

1.27 

- 

804 

1096 

/a 

1 

10 

5Rc 

•4.0 

1.76 

321.6 

485 

691 

h 

0 

a 

5  Re 

2.59 

1.58 

- 

604 

786 

fz 

6 

11 

•  5°F 

4.8 

1.75 

275.6 

449 

- 

fz 

- 

223 

20Rc 

2.54 

0.90 

612 

894 

- 

h 

3 

122" 

10R^ 

3.64 

1.27 

- 

662 

1006 

h 

9 

222 

20R5 

3.05 

0.90 

519 

830 

- 

h 

10 

12 

0^ 

5.9 

1.76 

231.8 

437 

730 

fz 

11 

10 

5^ 

4.0 

1.76 

321.6 

485 

691 

fz 

12 

130" 

10*7 

5.15 

1.27 

296 

644 

- 

fz 

13 

228 

20°? 

3.95 

0.90 

408 

828 

- 

h 

n 

122" 

10R5 

3.64 

1.27 

- 

622 

1006 

fz 

la 

12 

5R: 

5.9 

1.76 

231.3 

437 

730 

fz 

16 

227 

20R: 

4.60 

0.90 

353 

886 

- 

fz 

fz 

17 

226 

20°? 

5.15 

0.90 

309 

825 

- 

IS 

10 

4.0 

1.76 

321.6 

485 

691 

fl 

19 

11 

o'T 

4.8 

1 .75 

275.6 

449 

- 

fl 

20 

12 

5.9 

1.76 

231.8 

437 

730 

fx 

21 

130" 

lO^c 

5.15 

1.27 

296 

644 

- 

fx 

22 

227 

20  Re 

4.60 

0.90 

353 

886 

- 

fx 

23 

223 

20Rr 

2.54 

0.90 

612 

894 

- 

fx 

2  1 

228 

20re 

3.95 

0.90 

408 

828 

- 

fx 

2  a 

222 

20c7 

3.05 

0.90 

519 

830 

- 

fx 

20 

220" 

2  OR? 

5.15 

0.90 

309 

825 

- 

fx 

Measurements  by 

T  R.  Gururaja.  Materials  Research  Laboratory. 

The  Pennsylvania  State 

l  n  i  ve  r-  i  r  y . 

"Sample  also  tested  at  Stanford  University. 
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manufacturer’s  data,  using  piezoelastically-stiffened  elastic  constants  and  assuming  isotropy. 
Spurrs  epoxy  parameters  were  measured  at  Stanford,  using  pulse  echo  methods  to  determine 
the  velocities.  Frequencies  listed  in  the  table  are  the  anti-resonance  frequencies  for  the 
longitudinal  thickness  mode  (/J  and  the  resonance  frequencies  for  the  lateral  modes  (/2,/3 ). 
In  Fig  '-’(a),  the  resonance  frequencies  are  points  of  impedance  minimum  and  the  anti- 
resonance  frequencies  are  points  of  impedance  max'imum.  A  resonance  frequency  corresponds 
approximately  to  the  short-circuit  natural  frequency  of  the  freely  vibrating  piezoelectric 
resonator  mode,  and  an  anti-resonance  frequency  to  the  open-circuit  natural  frequency  of 
the  free  resonator  mode.  Since  our  composite  theory  does  not  yet  include  piezoelectricity,  the 
choice  of  which  experimental  frequency  to  compare  with  theory  is  not  immediately  obvious. 
Our  choice  was  based  on  the  following  considerations:  (1)  In  the  Mason  model  for  a  parallel 
plate  transducer/resonator,  the  open  circuit  (or  anti-resonance)  frequencies  correspond  to 
pure  n\/2  standing  waves  across  the  plate;8  for  this  reason,  the  anti-resonance  frequencies 
are  chosen  to  compare  with  the  theory  of  the  longitundinal  thickness  resonance  (/j  =  /(). 
j-J)  For  the  lateral  resonances  (/2,/3),  the  standing  waves  are  parallel  to  the  electrode  faces 
(Fig.  '2(b))  and  the  choice  between  resonance  and  anti-resonance  frequencies  is  not  clear  at 
this  point:  consequently,  the  resonance  frequencies  were  chosen  simply  because  they  give  a 
better  fit  to  the  theory. 

In  subsection  Q.A  above,  it  was  noted  that  the  constant  strain  model  of  Appendix  A  is 
more  appropriate  for  calculating  the  longitudinal  thickness  resonances,  while  the  constant 
stress  model  of  Appendix  B  is  more  appropriate  for  the  lateral  resonances.  Figures  6-11  com¬ 
pare  theory  and  experiment  for  these  two  cases.  Figures  6-8  represent  the  dispersion  curves 
for  the  symmetric  Lamb  modes  *  calculated  from  spatially  averaged  stiffnesses  (the  constant 
strain  model).  On  these  curves  the  frequencies  of  the  longitudinal  thickness  resonances  for 
samples  in  Table  I  are  plotted  at  3b  —  0.  They  lie  very  close  to  the  3b  =  0  intercept 
of  the  second  symmetric  Lamb  mode,  a  point  that  correponds  to  the  longitudinal  thickness 
resonance  in  a  homogeneous  plate.  Plots  of  lateral  resonance  frequencies  (at  3  =  2x/d,  cor¬ 
responding  to  the  resonance  conditions  in  Fig.  2(b)),  lie  far  below  the  Lamb  waves  spectrum 
in  Figs.  6-8.  This  is  unreasonable  physically  because  the  lateral  resonances  correspond  to 
Lamb  waves  that  are  Bragg  scattered  by  planes  of  piezoelectric  rods;  and  these  points  are 
therefore  not  shown. 

In  Figs.  9-11.  the  curves  are  Lamb  wave  spectra  calculated  from  spatially-averaged  com¬ 
pliances  ( i.e. ,  the  constant  stress  model).  The  experimental  longitudinal  thickness  resonances 
[fi )  are  now  seen  to  be  in  poor  agreement  with  the  theory,  but  the  lateral  resonances  (/2,  /3 ) 
(plotted  at  3  =  2 jr/d)  lie  close  to  the  lower  modes  of  the  symmetric  Lamb  spectrum,  as  ex¬ 
pected  from  Appendix  B.  The  experimental  points  are  distributed  along  the  dispersion-curves 
*  Antisymmetric  modes  cannot  be  excited  in  plates  driven  by  PZT  rods  poled  normal  to  the  plane  of  the  plate. 
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Lamb  Wave  Curves  for  Constant 
Strain  Average  Producers 

20%  volume  PZT-5A 
in  Spurrs  Epoxy 

=  3453.4  m/s 

V  >=  1890.4  m/s 
s 


20.0  JO  .0  40.0  SO.O 

3b 


FIGURE  8 

Same  as  Fig.  8  (20cc  volume  ratio) 


Comparison  of  constant  stress  model  with  experiment  (5^  volume  ratio). 
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and  trace  out  the  Lamb  wave  dispersion  curves  quite  accurately. 

C.  Discussion 

Experimental  points  plotted  in  Figs.  6-11  clearly  show  the  expected  behavior — i.e., 
longitudinal  thickness  resonances  are  best  represented  by  the  constant  strain  model  and 
lateral  resonances  by  the  constant  stress  model,  for  volume  percentages  of  PZT  as  large  as 
'2Ql'l .  In  this  comparison,  the  longitudinal  thickness  resonance  is  analysed  by  completely- 
neglecting  any  coupling  into  other  space  harmonics,  and  the  lateral  resonances  are  analysed 
by  the  coupled  mode  approximation.  These  approximations  are  expected  to  be  somewhat 
inaccurate,  even  in  the  .j rc  volume  ratio  case,  and  this  fact  may  account  for  some  of  the 
disagreement  observed  in  Figs.  6-11.  Continuing  analysis  should  add  more  and  more  space 
harmonics  to  the  calculation  until  there  is  no  further  change  in  the  result. 

Another  possible  source  of  error  is  the  incomplete  treatment  of  piezoelectricity  and 
anisotropy  in  the  PZT.  One  important  consequence  of  piezoelectric  coupling  was  discussed 
in  subsection  II. B.  with  respect  to  the  choice  of  the  longitudinal  thickness  anti-resonance 
for  comparison  with  theory.  The  theoretical  basis  for  choosing  to  compare  lateral  mode 
experiments  with  theory  on  the  basis  of  resonance  frequencies  has  not  yet  been  formulated. 
Another  area  of  uncertainty  is  how  best  to  characterize  piezoelectric  stiffening  of  the  elastic 
constants  in  a  composite  medium. 

Despite  these  unresolved  questions,  agreement,  of  theory  and  experiment  in  Figs.  6-11 
is  quite  sufficient  to  permit,  drawing  practical  conclusions  and  to  explain  and  predict  various 
aspects  of  the  behavior  of  composites.  In  many  experimental  resonance  curves  a  number 
of  small  supplementary  resonances  are  observed  above  and  below  fn  and  /(2  These  can 
be  explained  by  observing  that  Eq.  (3)  couples  among  all  Lamb  waves  in  all  orders  of  space 
harmonic,  and  that  additional  couplings  occur  in  the  two-dimensional  version  of  the  equation. 
Bragg  resonances  therefore  occur  for  the  higher  order  Lamb  modes,  but  are  weakly  excited 
by  the  electrodes. 

Understanding  the  resulting  structure  of  the  dispersion  curves  is  best  reached  by  using 
the  reduced  zone  scheme  for  the  dispersion  curves,  as  in  the  electronic  theory  of  solids.7  This 
can  be  visualized  by  referring  to  Fig.  4,  and  supposing  that 


which  corresponds  to  the  d  =  X  stopbands  discussed  above.  The  reduced  zone  scheme  is 
generated  by  folding  the  figure  along  the  vertical  line  through  dnb.  so  that  the  parts  of  the 
curves  on  the  right  cross  the  parts  of  the  curves  on  the  left.  Each  crossing  point  corresponds 
to  a  Bragg-matching  condition,  where  the  crossed  curves  split  into  stopbands.  Crossings 
on  the  vertical  line  through  3nb  =  2xb/d  represent  Bragg-scattering  of  a  Lamb  mode  into 
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itself;  these  are  the  most  strongly  excited.  Other  crossings  are  for  Bragg-scattering  from  one 
Lamb  mode  into  another.  At  these  points,  the  band  edge  resonances  do  not  have  piezoelectric 
rods  vibrating  in  phase,  so  that  they  are  less  strongly  excited  than  in  the  first  class  of  band 
edge  resonances.  Even  in  this  one-dimensionally  periodic  model,  there  exists  a  large  density 
of  stopbands  in  the  frequency  range  above  the  lowest  stopband  frequency;  and  in  a  two- 
dimensionally  periodic  structure  this  density  is  further  increased  by  the  proliferation  of  two- 
dimensional  Brillouin  zones.  Such  features  will  be  extremely  important  in  determining  the 
scattering  behavior  of  periodic  composite  materials,  and  may  even  be  applicable  in  designing 
and  fabricating  anti-reflection  layers. 

The  discussion  in  the  paragraph  above  is  also  relevant  to  the  behavior  of  observed 
resonance  spectra  curves  as  a  function  of  the  ratio  of  b/d  in  Fig.  4.  Since  the  curves  are 
normalized  to  6.  setting  3n  —  2 x/d  places  the  vertical  dashed  line  at  a  position  that  moves  to 
the  left  as  d  increases.  This  means  that  the  position  of  the  lowest  lateral  stopband  resonance 
moves  steadily  to  the  left  with  increasing  d,  as  can  be  tracked  by  comparing  the  points 
plotted  in  Figs.  0-11  with  their  listed  composite  dimensions  in  Table  I.  Eventually,  as  d 
increases,  the  lowest  stopband  frequency  approaches  coincidence  with  longitudinal  thickness 
resonant  frequency  (at  the  3nb  —  0  intercept  of  the  second  Lamb  mode  in  Figs.  6-8,  since  the 
constant  strain  model  must  be  used  for  this  resonance).  This  phenomenon  has  been  observed 
in  experiments  at  Penn  State,  and  an  apparent  coupling  of  the  thickness  and  lateral  modes  has 
been  found  to  occur  as  their  frequencies  cross.  Practically,  the  most  important  impact  of  this 
Lamb  mode  Floquet  theory  is  in  providing  a  physical  basis  for  optimally  choosing  the  ratio 
b/d  in  transducer  and  resonator  design.  However,  explaining  the  mode  coupling  phenomenon 
observed  at  Penn  State  may  also  turn  out  to  have  significant  practical  consequences. 
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ffl.  PMN-BASED  TRANSDUCERS 


The  goal  of  the  research  is  to  study  acoustic  wave  propagation  in  ceramic  materials 
in  which  spatially  varying  polarization  patterns  are  established,  with  the  aim  of  exploring 
novel  signal  processing  technology  and  applications.  Initial  work  will  be  conducted  with 
lead  magnesium  niobate  (PMN)  and  lanthanum  modified  lead-zirconate  lead-titanate  (PLZT) 
ceramics.  These  materials  exhibit  a  strain  which  is  proportional  to  the  square  of  an  externally 
applied  electric  field,  i.e.,  electrostriction.  In  addition,  suitable  compositions  do  not  possess 
remanent  polarizations.  It  is  this  latter  property  that  may  be  most  interesting  for  certain 
device  applications,  since  it  permits  rapid  switching  of  transducer  configurations. 

Three  experiments  are  currently  being  planned.  The  first  is  a  basic  study  of  some  electro¬ 
mechanical  properties  of  the  above-mentioned  ceramic  materials.  Since  these  materials  do 
not  possess  a  remanent  polarization,  the  electro-mechanical  parameters  are  a  function  of  an 
externally  applied  biasing  field.  This  first  experiment  will  measure  this  dependence  using 
a  plate  transducer  structure  of  thickness  t  and  area  VF  by  L.  The  parallel  plate  structure 
was  chosen  so  that  the  externally  applied  field  could  be  determined  from  readily  measured 
factors  (in  i h is  case  the  externally  applied  voltage  and  the  thickness  of  the  wafer).  The  basic 
electrical  measuring  circuit  is  shown  in  Fig.  12.  It  is  an  impedance  measuring  system. 

The  second  experiment  usts  the  structure  shown  in  Fig.  13(a).  It  is  intended  for  use  as 
a  surface  acoustic  wave  (SAW)  device.  The  idea  for  this  structure  is  based  upon  other  SAW 
devices  using  piezoelectric  materials  that  can  be  modeled  as  possessing  the  hexagonal  type 
of  symmetry  In  those  devices,  the  material  is  oriented  so  that  the  c-axis  is  perpendicular  to 
the  major  face  containing  the  surface  electrodes.  For  ceramics,  the  “c-axis”  will  be  colinear 
with  the  electric  field.  Thus,  in  our  case,  a  DC  bias  field  will  be  established  between  the 
full  surface  electrode  on  one  surface  and  the  electrode  pattern  on  the  opposite  surface.  The 
relationship  between  the  various  dimensional  parameters  is  governed  by  two  factors:  acoustics 
and  electrostatics.  Electrostatic  considerations  lead  to  the  following  relations: 

(w  +  g)/t  <  <  1  w l g  >  1  (5) 

Basically,  these  relations  insure  that  the  configuration  produces  an  electrostatic  field  which 
approaches  the  ideal  parallel  plate  field.  One  practical  consideration  is  that,  from  an  operating 
voltage  point  of  view,  we  do  not  want  t  to  become  too  large.  Acoustic  properties  lead  to  the 
following  relations: 

(u’  4-  g)/t  <<  I  w/g  =  1  (6) 

Clearly,  the  first  relation  is  desirable  since  we  want  pure  surface  waves,  which  require  that 
the  second  surface  be  very  far  from  the  top  surface. 

Finally,  we  need  to  study  the  effect  of  the  fringing  fields  near  the  electrode  h  eS  We 
intend  two  methods  for  this  study.  The  first  is  to  develop  a  computer  program  *  solve 


FIGURE  12 

PMN  impedance  measuring  circuit. 

the  electrostatic  field  problem,  and  the  second  is  to  experimentally  measure  the  field.  This 
second  method  will  use  the  PLZT  material.  This  ceramic  material  is  optically  active,  i.e.,  it 
is  uniaxially  birefringent.  Since  the  dielectric  properties  of  both  ceramics  are  very  similar, 
the  fields  obtained  for  the  PLZT  device  will  be  directly  applicable  to  the  PMN  device.  The 
structure  intended  for  this  experiment  is  shown  in  Fig.  13(b). 

We  plan  to  use  two  different  polarizer  configurations:  a  circular  polarizer  and  a  linear 
polarizer.  With  the  circular  polarizer  configuration  we  can  measure  the  magnitude  of  the 
eleerric  field  irrespective  of  its  direction.  Then,  using  the  linear  polarizer,  we  can  measure 
the  orientation  of  the  field.  Thus,  with  this  approach,  the  computer-generated  solution  can 
be  checked  and  confidence  built  for  its  use  in  other  configurations. 
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IV.  INTERFACE  WAVES 


Ferroelectric  ceramics  poled  with  a  pattern  of  spatially  varying  remanent  electric 
polarization  have  been  employed  in  a  number  of  novel  device  applications.  These  range 
from  optical  displays  and  memories  in  transparent  piezoceramics  to  Fresnel-type  focusing 
transducers  for  ultrasonic  imaging  in  medicine  and  nondestructive  evaluation.  For  these 
devices,  a  uniformly  poled  piece  of  ceramic  is  reverse-poled  in  certain  regions  to  effect  a 
desired  spatial  pattern  of  reversed  piezoelectric  coupling.  Transition  zones  separating  regions 
of  unlike  remanent  polarization,  analogous  to  domain  walls  in  ferroelectric  crystals,  arise  as  a 
consequence  of  this  regional  counter-poling  process.  While  the  transition  from  one  polariza¬ 
tion  state  to  another  in  a  single  crystal  is  quite  abrupt  (domain  wall  thicknesses  less  than  a 
lattice  constant  are  typical),  the  transition  in  a  differentially-poled  ceramic  is  comparatively 
gradual.  It  has  been  shown  theoretically  and  experimentally  that  both  abrupt  and  gradual 
polarization  transition  zones  are  capable  of  supporting  guided  acoustic  waves.8  Abrupt  tran¬ 
sitions  were  realized  by  glueing  together  two  uniformly  polarized  ceramic  plates  along  an  edge, 
with  the  two  polarizations  reversed.  Gradual  transition  zones  are  fabricated  by  depositing 
electrode  patterns  on  an  already  poled  ferroelectric  ceramic  and  counter-poling  under  the 
electrodes.8  This  creates  waveguiding  polarization  transition  zones  around  the  edges  of  the 
electrodes,  in  spatial  patterns  determined  by  the  electrode  geometries. 

In  this  way.  single  mode  acoustic  waveguides  can  be  fabricated  in  ferroelectric  ceramic 
wafers,  using  simple  inexpensive  differential  poling  techniques.  ‘‘Acoustostrip"  waveguide 
components  based  on  this  concept  may  be  conceived  in  all  of  the  microstrip  waveguide 
geometries  utilized  in  microwave  integrated  circuit  technology' — i.e.,  parallel  line  directional 
couplers,  ring  resonators,  hybrids,  channel-dropping  filters,  etc.  Directional  coupling  action 
has  been  already  demonstrated  experimentally — first,  at  1  \IHz9  and,  more  recently,  at 
several  megaherz.  The  significance  of  this  new  technology  is  that  it  has  potential  for  realizing 
integrated  acoustostrip  filters  and  signal  processing  systems  (operating  up  to.  perhaps.  10 
MHz)  on  a  single  ferroelectric  wafer.  A  possibility  also  exists  for  creating,  on  single  crystal 
wafers,  poled-in  systems  of  this  kind  for  operation  at  frequencies  up  to  hundreds  of  megaherz. 
M  aterials  such  as  PMN  (section  in,  above)  may  eventually  be  developed  in  sufficiently  high 
quality  to  permit  realization  of  switchable  acoustostrip  waveguide  systems  controlled  by  an 
electric  bias  field. 

A.  Waveguide  Fabrication  and  Characterization 

In  poling  an  acoustostrip  waveguide  duct  into  an  initially  uniformly  poled  ferroelectric 
ceramic  wafer,  a  variety  of  electrode  structures  may  be  used.  For  a  single  duct,  the  electrode 
covers  the  entire  half  of  the  wafer  on  one  side  of  the  duct;  for  a  pair  of  ducts,  a  strip 
electrode  is  used  (Fig.  11  (top)).  In  either  case  a  waveguide  duct  is  centered  on  each  zero 


electrode  electrode 
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FIGURE  14 


Aco'istostrip  waveguide  processing  on  a  ferroelectric  wafer 
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crossing  of  the  polarization  profile.  The  experimental  curve  at  the  top  right  of  Fig.  14  shows 
a  profile  creating  two  such  ducts,  spaced  approximately  by  the  width  of  the  strip  electrode. 
(Experimental  procedures  and  instrumentation  are  dicussed  in  section  IV. C  below).  Another 
electrode  configuration,  with  tapered  ends  suitable  for  directional  coupler  fabrication,  is 
shown  in  Fig.  15. 

Parameters  important  in  the  fabrication  of  these  waveguides  are:  (1)  properties  of  the 
ferroelectric  ceramic,  (2)  thickness  of  the  wafer,  and  (3)  curvature  of  the  electrode  edges.  The 
properties  of  the  ceramic  play  a  fundamental  role  in  establishing  the  width  of  the  waveguide 
duct,  which  is  related  to  the  slope  of  the  polarization  transition  profile  (Fig.  14,  top  right). 
This  width  is  important  because  it  determines  whether  or  not  the  waveguide  propagation  is 
single  mode  or  multimode.  Experiments  have  been  performed  on  both  a  hard  ceramic  (PZT-4) 
and  a  soft  ceramic  (PZT-5H),8  showing  that  PZT-4  produced  narrow  ducts  (singlemode)  and 
PZT-oH  broad  ducts  (multimode). 

A  conflicting  requirement  on  the  choice  of  ceramic  is  that  it  must  have  as  large  a  Q  as 
possible  at  megaherz  frequencies.  There  exists,  at  present,  essentially  no  qualitative  data  in 
the  ferroelectric  ceramics  literature  concerning  Q  as  a  function  of  frequency.  However,  our 
waveguide  experiments  show  that  PZT-5H  (in  which  single  mode  waveguides  were  easily  real¬ 
ized)  has  unacceptably  large  attenuation  losses.  A  compromise  candidate  has  been  found  in 
PZT-8.  This  is  a  hard  ceramic  with  good  high  frequency  Q,  but  it  also  generates  counter-poled 
ducts  of  widths  only  slightly  broader  than  those  in  PZT-oH.  The  exact  physical  phenomenon 
determining  the  duct  width. as  not  yet  been  identified;  but  in  the  materials  listed  a  correlat  u 
has  been  noted  between  large  elastic  stiffness  and  broad  duct  width.  Barium  titanite  ceramics 
should  also  be  studied  in  the  future,  to  further  explore  this  question. 

Finally.  the  thickness  of  the  ceramic  wafer  and  the  curvature  of  the  electrode  edge 
influence  the  width  of  the  electrode  fringing  field,  and  therefore  also  control  the  width  of  the 
region  subject  to  depoling.  (Figure  16  illustrates  a  calculated  plot  of  the  fringing  field  near 
the  edge  of  a  depoling  electrode  on  ceramic).  The  width  of  the  fringing  field  region  decreases 
with  the  thickness  of  the  wafer;  and  it  has  been  confirmed  experimentally  that  duct  widths 
are  narrower  on  thinner  wafers.  This  is  another  factor  to  be  considered  in  designing  these 
acoustostrip  waveguide  structures. 

Measurements  of  guided  wave  propagation  were  reported  on  both  PZT-4  and  PZT-5H 
single  acoustostrip  waveguides  in  Reference  8.  These  experiments  were  performed  using  a  new 
method  of  wave  excitation  (edge  bonded  transducers)  that  has  turned  out  to  be  a  superior 
method  for  creating  clean,  localized  acoustic  excitation  profiles.  Laser  probe  measurements 
were  made  in  the  pulse  mode  to  determine  the  amplitude  profile  of  the  propagating  duct 
wave.  Very  accurate  propagation  velocity  measurements  were  also  carried  out,  and  these  gave 
glue-bond  and  transition  wave  velocities  in  very  good  agreement  with  theoretical  predictions. 
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FIGURE  15 

Differential  poling  apparatus 
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FIGURE  18 

Calculation  of  fringing  fields  outside  poling  electrodes. 
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Improper  poling  conditions  and/or  inferior  ceramic  may  result  in  cracking  within  the 
poling  transition  zone.  Interface  waves  scatter  strongly  from  such  flaws  and  provide  an 
excellent  means  for  probing  the  integrity  of  transition  zone  structures  intended  for  waveguide 
or  other  applications. 

B.  Directional  Couplers 

The  bottom  part  of  Fig.  14  shows  our  first  acoustostrip  waveguide  structure  designed  to 
test  the  directional  coupling  principle.9  Two  closely-spaced  waveguide  ducts  were  created  in 
the  depolarized  transition  zones  at  the  edges  of  a  strip  electrode  (Fig.  14,  top).  Coupling 
between  the  waveguides  is  adjusted  by  controlling  the  electrode  width,  so  that  the  two 
guided  wave  profiles  overlap.  In  the  light  of  the  discussion  in  section  IV.A  above,  it  can 
be  seen  that  this  overlap  will  also  be  affected  by  the  properties  of  the  ceramic  and  by  the 
wafer  thickness.  Directional  coupling  action  was  demonstrated  by  laser  probing  the  acoustic 
vibration  amplitudes  on  the  two  parallel  acoustostrip  guides.  An  example  of  measured  data 
is  shown  at  the  lower  right  of  Fig.  14,  where  the  wave  is  launched  at  the  top  of  the  set  of 
curves  and  proceeds  downward.  It  is  seen  that  transfer  to  the  second  guide  occurs  after  a 
travel  distance  of  26.6  mm.  This  data  was  taken  using  8-cycle  tonebursts  at  1  MHz  on  a 
PZT-5H  plate  3.81  mm  thick  and  with  a  diameter  of  approximately  51  mm.  As  described 
in  Reference  9,  the  effective  spacing  of  the  coupled  waveguides  was  subsequently  decreased 
by  applying  a  reversed  depoling  voltage  (6.3  KV/cm)  to  the  electrode  for  60  seconds.  The 
effect  of  this  was  to  increase  the  coupling  and  to  decrease  the  energy  transfer  distance  from 
26.6  mm  to  approximately  11.6  mm.  Voltage  control  of  directional  coupling  could  be  very 
significant  in  PMN  structures,  if  waveguide  fabrication  should  prove  to  be  practical  in  these 
materials. 

A  new  generation  of  acoustostrip  directional  couplers  has  now  been  fabricated,  with 
tapered  duct  spacing  at  the  input  and  output  to  provide  better  isolation.  Figure  15  shows 
one  of  these  devices  in  the  poling  station,  where  the  tapered  depoling  electrode  is  clearly 
visible.  In  Fig.  17,  the  processed  device  is  shown  with  the  poling  electrode  still  in  place; 
and  in  Fig.  18  it  is  shown  mounted  in  the  text  fixture,  with  the  electrode  removed.  These 
devices  are  being  tested  in  the  pulse  mode  at  frequencies  up  to  7  MHz.  The  objectives  are  to 
carefully  quantify  the  S-parameters  of  the  devices  and  to  obtain  quantitative  data  relating 
performance  to  the  fabrication  parameters  (i.e.,  poling  profile  shape,  wafer  thickness,  and 
material  properties).  This  will  lead  to  better  physical  understanding  of  how  the  polarization 
response  of  the  ceramic  can  be  tailored  to  suit  various  eventual  device  applications  (directional 
couplers,  resonators,  filters,  etc  ). 

C.  Instrumentation 

To  obtain  the  above  results  and  to  prepare  for  achieving  the  device  characterization  goals 
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FIGURE  18 

PZT-5H  directional  coupler  in  test  fixture.  Poling  electrodes  have  been  removed 
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stated  at  the  end  of  section  IV. B,  it  has  been  found  necessary  to  improve  our  laser  probe  and 
poling  system. 

Laser  Probe.  During  the  current  reporting  period,  the  laser  probe  was  interfaced  with 
a  dedicated  laboratory  microcomputer.  Software  was  written  and  hardware  developed  to 
control  the  sample  scanning  mechanisms  fully  automatically,  and  to  acquire  and  plot  data. 
(This  facility  has  also  played  an  essential  part  in  the  composite  material  studies  of  section  I) 
The  electronics  of  the  laser  probe  detection  system  has  been  modified  to  permit  measurement 
of  acoustic  tone  bursts  at  carrier  frequencies  up  to  7  MHz. 

Poling  Apparatus.  To  permit  rapid  characterization  of  the  poling  conditions,  new  record¬ 
ing  inst  rumentation  has  been  developed.  The  circuit  used  is  illustrated  in  the  insert  of  Fig.  19. 
where  the  series  resistor  is  intended  to  limit  current  flow  in  the  event  of  breakdown,  and  the 
current  integrator  monitors  the  degree  of  poling  achieved  under  various  conditions.  During 
a  poling  operation,  the  amount  of  charge  flowing  through  the  integrator  is 

•^Q  total  =  initial  charge  T  AQ  discharge 

I 

=  A(AP), 

where  A  is  the  electrode  area  and  AP  is  the  change  in  ceramic  polarization  due  to  poling. 
Figure  19  illustrates  data  taken  in  this  way  for  a  PZT-5H  sample,  with  the  process  allowed 
to  run  to  completion  at  each  value  of  poling  field.  The  same  system  can  be  used  to  quantify 
poling  conditions  as  a  function  of  time. 

V.  SUMMARY  OF  SIGNIFICANT  ACOMPLISHMENTS 

(1)  Floquet  theory  of  periodic  composite  plates  developed  and  compared  satisfactorily 
with  frequency  resonance  spectra  measured  at  Penn  State  on  a  large  variety  of 
samples.  This  theory  predicts  dependence  of  the  lateral  resonance  frequencies  on  the 
thickness/spacing  ratio  for  the  PZT  rods  in  plastic  plates,  with  percentage  volume 
rates  up  to  20cc. 

(2)  Characterized  the  dependence  of  interface  waveguide  duct  thickness  (acoustostrip 
waveguide)  on  the  material  properties  of  the  ceramic  and  on  the  plate  thickness. 

(3)  Experimental  demonstration  of  acoustostrip  directional  coupler  action  at  frequencies 
up  to  7  MHz. 

(-1)  Extended  the  operation  of  the  laser  probe  up  to  7  MHz  in  the  pulsed  mode. 

(5)  Developed  new  circuitry  for  rapidly  monitoring  the  poling  state  in  fabricating 
differentially  poled  samples. 


32 


VI.  REFERENCES 


1  T.R.  Gururaja.  W.L.  Schulze,  L.E.  Cross,  B.A.  Auld,  Y.A.  Shui,  and  Y.  Wang,  “Resonant 
Modes  of  Vibration  in  Piezoelectric  Composites  with  Two  Dimensional  Periodicity," 
presented  at  the  5th  European  Meeting  on  Ferroelectricitv,  Benalmadena  (Malaga), 
Spain,  September  1983  (Ferroelectrics  54,  pp.  183-186  (1983)). 

2.  B.A.  Auld,  Y.  Shui,  and  Y.  Wang,  “Elastic  Wave  Propagation  in  Three-Dimensional 
Periodic  Composite  Materials,”  presented  at  the  International  Conference  on  the 
Dynamics  of  Interfaces,  Lille,  France,  September  1983  (to  be  published  in  Journal  de 
Physique). 

3.  B.A.  Auld.  H.A.  Kunkel,  Y.A.  Shui,  and  Y.  Wang,  “Dynamic  Behavior  of  Periodic 
Piezoelectric  Composites,”  Proe.  1983  IEEE  Ultrasonics  Symposium  (Catalog  No. 
83CH1947).  55-4-558  (1984). 

4.  S.  Ikegami.  I.  I'eda  and  S.  Kobavashi,  “Frequency  Spectra  of  Resonant  Vibration  in  Disk 
Plates  of  PbTiOs  Piezoelectric  Ceramics,”  J.  Acous.  Soe.  Am.  55.  339-344  (1974). 

5.  B.A.  Auld,  Acoustic  Fields  and  Waves  in  Solids,  Vol.  II,  Chapter  11  (New  York:  Wiley- 
Interscience,  1973). 

6.  B.A.  Auld,  Acoustic  Fields  and  Waves  in  Solids,  Vol.  I,  p.  324-327  (New  York.  Wiley- 
Interscience,  1984). 

7  J.M.  Ziman,  Principles  of  the  Theory  of  Solids,  p.  23  and  pp.  72-73  (Cambridge: 
University  Press,  1964). 

8.  H.A.  Kunkel  and  B.A.  Auld.  “The  Use  of  Guided  Acoustic  Interface  Waves  in  the  Study 
of  Differentially  Poled  Piezoelectric  Ceramics,”  presented  at  the  1983  EEEE  International 
Symposium  on  APPLICATIONS  OF  FERROELECTRICS  (ISAF).  Gaithersburg,  Md„ 
May  1983  ( Ferroelectrics  51,  99-104  (1983)). 

9  H.A.  Kunkel  and  B.A.  Auld,  “Acoustic  Interface  Waveguide  Directional  Coupler,"  Proc. 
1983  IEEE  Ultrasonics  Symposium  (Catalog  No.  83CH1947),  562-565  (1984). 


34 


VH.  APPENDIX  A 


From  Proceedings  of  l  Itrasonics 
Symposium.  1983  (Reference  3) 

DYNAMIC  3EHAVI0R  OF  PERIODIC  PIEZOELECTRIC  COMPOSITES 


3.  A.  Auld,  H.  A.  Kunxel,  Y.  A.  Shui,  and  Y.  Wan g 


Eduard  L.  Ginzcon  Laboratory 
Scanford  University 
Stanford,  California  94305 


Abstract 

Comp  sice  piezoelectric  transducer  materials 
consisting  of  two-  and  three-dimensional  periodic 
arrays  of  piezoelectric  ceramic  elements  (i.e., 
rods,  soheres.  etc.)  in  a  plascic  macrix  now  pLav 
jn  tmoortanc  role  in  hydrophone  technology.  These 
materials  uere  conceived  and  developed  on  che  basis 
oc  scriccly  stacic  elasticity  considerations,  suit¬ 
able  cor  general  aspects  of  hydrophone  performance, 
puc  in  order  to  completely  understand  and  exploit 
such  materials  it  is  essencial  co  fully  develop  a 
theorv  of  their  dynamic  behavior.  This  paper  pre¬ 
sence  a  general  Floquet  theory  of  uave  propagation 
in  periodic  elastic  composites  and  develops  de¬ 
tailed  results  for  composite  places  consisting  of 
..erumic  rods  imbedded  in  a  plascic  macrix.  Laser 
prooe  measurements  of  the  vibracion  discribucion 
.a  such  places  show  good  agreement  vich  cheory. 

1.  Introduction 

The  design  and  fabrication  of  artificial  mate¬ 
rials  for  special  purposes  is  a  rapidly  growing 
tecnnology,  wich  applications  spanning  che  range 
-rom  mechanical  structures  to  electronic  devices. 

An  important  class  of  these  new  materials  is  che 
family  of  piezoelectric-polymer  composite  trans¬ 
ducers  for  various  sonic  and  ultrasonic  devices.  - 
In  designing  such  macariais,  consideration  is  given 
to  choosing  che  component  phases  of  the  composite 
material  so  as  co  maximize  some  figure  of  meric  for 
t.me  application  envisaged.  For  example,  in  an  elec¬ 
tromechanical  transducer  one  might  wane  co  maximize 
tne  piezoelectric  coupling,  minimize  the  dansicy  and 
optimize  the  acoustic  match  to  wacer.  Single  phase 
-acerials,  such  as  ceramic  piezoeleccrics  and  plas¬ 
tic  p  lezoeieccrics,  are  not  themselves  optimum  in 
a*,  these  respects  and  che  goal  In  composite  mate¬ 
rial  aesign  is  co  find  the  Pest  combination  of 
onases.  For  examole.  it  has  been  found  chat  com- 
’osites  consisting  of  ?ZT  rods  imoedded  in  a  plastic 
macrix  have  certain  piezoelectric  coupling  coer'fi- 
-.ents  chac  are  much  suoerior  to  those  of  homoge¬ 
neous  ?ZT. 1  Many  different  regular  two-phase  \PZT 
inu  Uastic)  composite  structures  (or  laccices)  nave 
teen  fabricated  and  scudied  in  che  low  frequency 
long  wavelength)  regime,  wnere  the  periodicity  of 
me  structure  is  small  tomoared  with  che  acoustic 
oavelengcn.  This  is  appropriate  co  hydrophone 


applications.  More  recently  interest  has  been 
directed  coward  che  behavior  of  these  materials  in 
che  frequency  range  where  the  structural  period  is 
comparable  to  a  wavelength  and  where  the  material 
itself  exhibits  frequency-dependent  behavior  simi¬ 
lar  co  chat  of  vibrations  of  a  crystal  lactice  near 
che  Brillouin  zone  boundary.  Quantitative  under¬ 
standing  of  this  phenomenon  is  necessary  for  an 
evaluation  of  che  role  of  chese  materials  in  trans¬ 
ducer  applications  at  higher  frequencies. 

2.  Floquec  Theory  for  Two-Dimensional  Composites 

Although  che  theory  presented  here  is  appli¬ 
cable  co  full  three-dimensional  lattices  of  arbi¬ 
trary  inclusions  imoedded  in  a  .macrix  and  can  be 
extended  co  periodic  composites  wich  three  or  more 
phases,  che  discussion  will  be  restricted  to  two- 
dimensional  structures  such  as  those  shown  in 
Fig.  1.  The  creacmenc  is  completely  parallel  co 
the  cheory  of  che  band  structure  of  solids-  but  it 
involves  vector,  rocher  chan  scalar,  wave  functions. 
Another  analogy  may  be  drawn  wich  che  vibracions  of 
a  crystal  lattice,  with  che  difference  that  the  com¬ 
posite  is  a  continuous  medium  racher  than  a  lactice 
of  point  masses  coupled  by  springs. 


SI  SI 


Figure  —  Two-dimensional  composite  unit 
cell  geometries. 

In  a  homogeneous  elastic  medium,  where  the 
elastic  constants  and  density  are  apaciail -  con¬ 
stant,  wave  propagation  at  frequency  -  -s  gov¬ 
erned  pv 
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vnere  tne  j-ae rivacive  actually  appears  tc  the  Left 
nr  che  elastic  conscanc  but  can  oe  shifted  to  the 
rignc-nana  side  wnen  the  constanc  is  uniform  in 
space.  This  step  cannoc  be  caicen  in  a  composite, 
where  spatiai  derivatives  of  the  elascic  constants 
oust  oe  considered.  To  do  this,  it  is  convenient 
to  write  Eq  .  . i)  as 


, . .  u  . 

wj  J 


C-) 


jnere  the  operators  are  defined  in  Ref.  3,  and  to 
separate  the  -material  parameters  into  spatial  aver¬ 
age  and  spacially  varying  parts. 


1 (x.y) 


CIJ  + 

o  +•  lo(x,yj 


(3a) 

(3b) 


suoscitution  reduces  Eq .  (2)  co 


•  x(3  -  Em'i/di 

—  x 

-  y(2  *  2-a/d) 


C  7b  > 


"or  conciseness,  summation  signs  are  omitted  ana 
summation  over  repeated  subscribes  is  understood. 
The  first  two  suosenpes  in  Eq  •  "a)  are  integers 

“  to  m-  *>)  describing  the  order  of  the  soace 
harmonid  and  the  third  subscript  i  is  a  vector 
polariaation  laoel.  Since  each  space  harmonic 
, ia)  is  a  three-dimensional  vector  field  it  must 
be  described  by  three  vector  components,  these 
could  be  taxen  as  components  with  respect  to  the 
coordinate  axes.  However,  it  is  much  more  useful 
co  take  as  a  basis  the  polarization  vectors  of 
elastic  waves  propagacing  in  che  average  medium 


with  wave  vector 


Urn 


Only  isocropio  media  will 


be  considered  here  and,  in  chis  case,  the  decom¬ 
position  is  as  shown  w  rig.  2.  this  means  that 


__  _  7 

(.  -  biw"t)  •  u  -  P(x,y)  •  u 


i) 


a  3  «  3  macrix  equation,  the  first  term  on  the 
-eft  is  the  operacor  on  che  left-hand  side  of  Eq . 

21 ,  written  for  the  spacial  average  elascic  con¬ 
stants,  ana  the  second  cerm  is  proportional  co  che 
teen city  macrix.  On  the  rlghc. 


? ( x , y )  »  -  (If  +  leu  I )  - 

a  '  a  a  a 


(5) 


wnere  the  quantity  in  parencheses  has  the  average 
conscancs  appearing  in  che  left-hand  side  of  Eq . 

i)  replaced  by  the  spatially  varying  values  in 
Eq  !3).  The  second  term  cakes  account  of  the 
spacial  derivatives  of  the  elastic  conscancs.  It 
las  the  form 
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Figure  2  —  Vector  soace  harmonica  of  order  lm  for 
a  periodic  elaacic  composite.  <,a)  Longi¬ 
tudinal  polar izacion .  (b)  Vertical  shear 

polarization .  <  c)  Horizontal  snear 

polarization . 


the  ooianzacion  vectors  in  Eq .  ai  are  eigenvec¬ 
tors  of  the  left  side  of  Eq .  - 1  .  3ecau  se  of  :ms, 
Ed.  .**1  reduces  to  an  infinite  set  of  touoied  ecua- 
tions  for  tn«  soace  harmonic  aaniitudea 


.  pc: 


wnere  tne  : 
e.as  t ic  wave 


:n  the  . elt-nanc  side  are  average 
'.starts  The  :  o  u  o  *  i  r.  a  :  d  n  s  c  a  n  t  i  are 


tr.e  square  rad  laccice  of  Fig.  1(a).  For  che  Cir¬ 
cular  rad  geomecrv  the  smc  functions  are  replaced 
cv  ;  inc  functions.  A  dispersion  reiacion  between 
and  3  in  Fq.  (7b)  13  obcained  by  setting 
::e  decermmanc  of  Fq  .  (3)  ;o  zero,  performed  in 
practice  ov  cruncacir.g  the  system  and  adding  addi¬ 
tional  terms  until  satisfactory  accuracy  is  at¬ 
tained.  In  the  limit  of  D  »  0  in  Fig.  I  the 
coupling  constants  go  to  zero,  and  the  solutions 
oecome  the  three  plane  elastic  solutions  traveling 
m  the  direction  of  3^  in  a  homogeneous  medium. 
Consequently,  the  solutions  (elascic  31och  vaves) 
of  Fq.  (3)  are  of  three  distinct  types,  consiscanc 
vi th  this  limiting  behavior. 


3.  Wave  Propagation  in  Special  Directions 

'Jnderscanding  the  wave  behavior  described  bv 
the  theory  above  and  applying  it  co  Che  interpre¬ 
tation  of  experiments  is  facilitated  by  considering 
special  propagation  directions.  If  3Z  is  zero 
.n  Fq .  1 7b),  the  assumed  solutions  are  z-indeoen- 
uenc  and  the  macrices  in  Eq .  (4)  partition  into  a 
3  <  2  matrix  involving  ux,uy  and  a  scalar  part 
involving  uz  .  This  is  easily  pictured  physi¬ 
cally  if  the  generation  of  space  harmonics  is 
visualized  as  arising  from  multiple  scattering  of 
xv-oropagacing  vaves  bv  a  periodic  square  arrav  of 
the  z-onenced  rods  in  Fig.  1;  and  It  means  chac 
z-oolarized  'and  xv-polarized)  solutions  can  be 
treated  independently. 

The  z-polarized  solution,  being  scalar,  has 
the  same  fora  as  an  electronic  wave  function  and 
is  therefore  particularly  useful  for  picturing  the 
pnvsical  benavior  of  these  vaves.  A  basic  feature 
.s  the  existence  of  stODbands,  vhich  can  be  related 
to  3ragg  scattering  from  planes  of  rods  having 
various  orientacions  in  the  square  lattice.  The 
essential  pnenomena  can  be  explained  by  reference 
to  Fig.  3.  Suppose  In  Fig.  3(a)  that  a  z-polarized 
<-oropagating  snear  vave  is  incident  on  the  grating. 
At  a  frequency  for  vhich  the  rod  spacing  d  is 
one-half  shear  wavelength,  constructive  reflection 
3ragg  scattering)  occurs  from  the  vertical  planes 
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Figure  3  -  time  cell  georaecries  and  laser  scan 
pachs.  (a)  Standing  vaves  at  the 
second  scopoand,  along  che  unit  ceil 
edge.  b)  Scanding  vaves  at  the  second 
stoooand,  along  the  anic  cell  diagonal. 


of  roas,  and  resonance  occurs  beeveen  adjacent 
vertical  planes.  The  x-prooagacing  vave  ts  also 
scattered  from  planes  of  rods  at  -3  oegrees  to  the 
x  and  y  axes;  but,  in  tnts  case.  3ragg  conci- 
tions  are  not  satisfied  and  resonance  toes  not 
occur.  A  scoooand  exists  sc  o  -  -5,1  tor  an  x- 
prooagating  vave,  but  tnere  is  no  resonanc  scat¬ 
tering  into  /-propagating  vaves.  It  t.ne  scoppanb 
edge  the  solution  is  s  scanding  vave  along  x  . 

If  the  scattering  is  veax,  the  vave  form  is  ac¬ 
curately  described  by  oniy  evo  space  na monies , 
those  for  vhich  the  corresoonding  left-nano  sices 
of  Eq .  (8)  are  simultaneously  near  resonance 
(0  0  and  -l  0;  vith  3  i  3X  )  .  From  t.ne  sauare 
symmecry  of  the  iaccice  a  similar  scanding  vave 
solution  along  y  ,  aesenbed  spproximateiv  o v 
space  harmonics  '0  3  and  3  -i;  vith  3  =  :  •, 

also  exists.  This  behavior  corresoonds  to  the  first 
stopband  of  the  lattice. 

When  the  rod  spacing  d  is  one  shear  vave- 
lengch,  a  comolacelv  different  scanding  vave  pat¬ 
tern  results,  as  illustrated  in  Fig.  3(S/.  In  this 
case  a  z-polarized  x-prooagacing  vave  is  again  3ragg 
scattered,  vith  a  phase  shift  of  I”  from  one  ver¬ 
tical  plane  of  rod3  to  an  adjacent  plane.  However, 
ic  is  easily  verified  chac  che  vave  is  also  3ragg 
scattered,  at  the  same  frequency,  by  the  -5°  Planes 
of  rods,  so  thac  resonanc  scattering  of  the  x-orooa- 
gating  vave  into  a  v-propagatmg  vave  occurs.  These 
y-propagacing  vaves  also  experience  resonance  scat¬ 
tering  of  the  same  cvpe.  As  a  result,  the  solution 
is  a  cvo-dimensional  resonant  standing  vave  pactem 
described,  for  veax  scattering,  oy  tour  space  har¬ 
monics  (0  0,  -2  0,  -l  l,  -l  -l,  vith  ;t  I  ix  ). 
Different  comoinaeions  occur,  consistent  vith  the 
symmecry  of  che  lattice.  Figure  3(a)  snows  the  one 
for  which  the  rods  sll  vibrate  in  Phase.  This  and 
che  ocher  two-dimensional  standing  wave  oacterns 
correspond  to  che  second  scopoand  of  the  square 
composite  lattice. 

Similar  comoencs  can  be  made  for  xv-polarized 
vave  propagation,  starting  from  an  x-ooiarized  x- 
propagating  or  a  y-poiarized  x-oropag3tin.g  vave. 

In  this  case  the  3ragg  scattering  argument  snows 
thac  only  longitudinal  soace  harmonics  or  s-esr 
space  harmonics  C3n  be  simultaneous!.'  -.ear  reso¬ 
nance,  a  conclusion  that  is  easily  confirmed  :v 
Fq.  (3). 

Anocher  special  case  is  the  one  :or  .-men 


x  / 

in  Fq .  v'b).  Suopose,  for  examo.e.  a  r-oc.ar.cec 
z-propagatmg  vave  -s  incident  on  the  .attics  ' 
3ragg  scattering  can  occur  aecause  :-e  structure  is 
uniform  along  z  However,  a  spates.  '.ci..a:icr 
of  che  pnase  fronts  is  produced  ere  utterance 
in  the  prooerties  of  Che  rocs  and  ere  ra c t  ..< .  Icr.- 
sequenciy.  soace  harmonics  are  gener3tec,  out  -ere 
are  resonant.  3i.mi.ar  conclusions  are  oocair.ee  :;r 
z-?ropagacmg  vaves  of  other  ?o  iar  icat  icr.s .  arc  ... 
are  consistent  v-.cn  the  system  or  oouo.ee  ,oacs  vr- 
momc  equations  Fq.  3V. 
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VUI.  APPENDIX  B 


LAMB  WAVE  FLOQUET  THEORY  FOR  A 
PERIODICALLY  COMPOSITE  PLATE 

Floquet  theory  for  an  infinite  isotropic  composite  medium  was  formulated  in  Appendix 
A  by  expressing  solutions  to  the  acoustic  wave  equation  in  a  series  of  plane  wave  harmonics 
Each  order  of  space  harmonic  (for  example,  the  fm-th  space  harmonic)  includes  plane  waves 
with  the  three  possible  plane  wave  polarizations  corresponding  to  the  space  harmonic  propaga¬ 
tion  factor  3  ^  (Fig.  2  of  Appendix  A).  There  are  therefore  three  space  harmonic  amplitudes 
aim,)  in  =  L.S.S')  in  each  order  (fm)  of  space  harmonic.  Substitution  of  the  Floquet  solu¬ 
tion  (or  space  harmonic  expansion)  into  Eq.  (2)  of  Appendix  A,  and  use  of  the  orthogonality 
properties  of  the  space  harmonic  wave  functions,  lends  to  an  infinite  set  of  coupled  equations 
for  the  space  harmonic  amplitudes  (Eq.  (8)  in  Appendix  A).  Standard  methods  can  be  used 
to  solve  these  equations,  including  the  coupled-wave  approximation  (in  which  only  space  har¬ 
monics  that  are  near  "resonance"— i.e. ,  have  the  right-hand  member  of  Eq.  (8)  near  zero — are 
retained).  Difficulties  arise  in  treating  the  composite  plate  problem  by  this  method,  because 
space  harmonics  are  coupled  both  by  the  spatial  inhomogeneities  (Ap,  Ac/M  and  by  the  plate 
boundary  conditions.  These  boundary  conditions  couple  not  only  space  harmonics  of  different 
orders  but  also  different  polarizations  in  the  same  order  of  space  harmonic.  Selection  of  the 
important  terms  in  the  Floquet  solution  therefore  becomes  very  complicated. 

The  difficulties  stated  above  can  be  avoided  by  taking  the  space  harmonic  terms  ..  the 
Floquet  solution  for  the  plate  to  be  Lamb  mode  and  SH  mode  wave  functions,  rather  than 
combining  infinite  plane  wave  space  harmonics  to  match  the  plate  boundary  conditions.  By 
choosing  Lamb  and.  SH  mode  wave  functions  as  a  basis,  one  pre-matches  the  space  harmonics 
to  the  plate  boundary  conditions.  To  illustrate  the  method,  a  simplest  example  (the  one- 
dimensionally  periodic  composite  plate  of  Fig.  4)  has  been  chosen.  For  plate  mode  Floquet 
analysis  the  acoustic  field  equations 
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arp  used,  rather  than  Eq.  fL)  of  Appendix  A.  Here,  the  spatial  average  dens.<y  p  and  com¬ 
pliance  a  are  defined  in  a  similar  wav  to  p  and  c  of  Appendix  A.  The  general  form  of  the 
Floquet  solution  is 
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where  u(r).  71  r)  are  functions  of  period  d  (expanded  in  a  Fourier  series),  and  an  are  the 
space  harmonic  amplitudes.  At  the  stress-free  surfaces  of  the  plate  =  0,6)  the  boundary 
conditions  are 

b  T  =y"  anb  ■  Tn  exp(—iJnx) 

n 

Since  this  condition  applies  for  all  values  of  x, 

b-T„=0  (3) 

for  all  n. 

As  in  Appendix  A,  the  individual  space  harmonic  wave  functions  on  the  right-hand  side 
of  Eq.  (2)  must  be  expanded  in  suitable  basis  functions.  In  a  structure  with  one-dimensional 
periodicity  and  propagation  parallel  to  the  periodicity,  the  Lamb  and  SH  modes  are  decoupled 
by  symmetry.  The  Floquet  solutions  may  therefore  be  taken  to  be  of  either  Lamb  mode  or 
SH  mode  species.  Here,  the  first  case  is  chosen;  and 


where  vnL,.  Tnl,  is  the  Lamb  mode  function  that  satisfies  the  acoustic  field  equations  at  and 
~  ’ ru,  (Fig.  4).  These  mode  functions  can  be  shown  to  satisfy  the  orthogonality  condition 

of  Eq.  (11.7-5)  in  Reference  5.  To  develop  coupled  space  harmonic  amplitude  equations,  Eq. 
(4)  and  Eq.  (2)  are  substituted  into  Eq.  (1).  This  result  is  reduced  to  a  set  of  coupled  space 
harmonic  equations  by  following  the  procedure  used  in  deriving  Eq.  (11.71)  of  Reference  ■> 
(taking  solution  "1"  to  be  Eq.  (2)  above  and  solution  "2”  to  be  the  Lamb  mode  wave  function 
with  n  =  rn  and  u  —  /<  (^'  =  J  in  Eq.  (4).  The  result,  corresponding  to  Eq.  (11.71), 
is  then  integrated  through  the  plate  thickness.  Expanding  the  divergence  operator  on  the 
right-hand  side  of  Eq.  (11.71),  applying  the  boundary  conditions  of  Eq.  (3),  using  the  fact 
that  the  rn.fi  solution  satisfies  the  field  equations  at  -c  =  and  integrating  over  the  unit 

cell,  finally  yields  the  coupled  space  harmonic  equations 

(^  )u  m*j  ==  ^  R  m/i.nt' ®  ne  ■  ( 5 ) 

riv 

The  coupled  constants  on  the  right-hand  side  of  Eq.  (5)  are  obtained  by  integrating  the 
A p.  As  matrix  in  Eq.  (1)  (weighted  by  multiplying  with  appropriate  combinations  of  Lamb 
mode  space  harmonic  functions)  over  the  unit  cell  of  the  composite. 
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